ytestgmmsvm_error = zeros(48, 1);
ydevgmmsvm_error = zeros(48, 1);

for k = 2 : 48
    ydevgmmsvm_error(k) = calc_error(gmmsvm_classify(Xdev, gmmdata2, HP, k), ydev);
    ytestgmmsvm_error(k) = calc_error(gmmsvm_classify(Xtest1, gmmdata2, HP, k), ytest);
end